On the Rate of Spontaneous Magnetic Reconnection 
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Magnetic reconnection is a topological rearrangement of the magnetic field lines, leading to the 
release of magnetic energy, which is thought to be associated with solar flares, coronal mass ejections 
and magnetospheric storms [l|, 01 • Despite magnetic field lines are supposed to be frozen into the 
well-conducting plasma, the reconnection observed in nature is, typically, fast, so that the rate of 
convergence of the magnetic field lines is the fraction of the Alfven speed, va- The Sweet-Parker 
solution P, 0] for the laminar thin current sheets predicts reconnection rates which are negligible for 
the solar or astrophysical conditions, this have prompted research into coUisionless reconnection 0- 
0] . The stochasticity of magnetic field lines due to ambient turbulence leads to fast reconnection 
and the rate was predicted to be proportional to kinetic energy density of ambient turbulence. Also, 
tearing instability of the thin current sheet [lol.[T]| was proposed as a driver of resistivity-independent 
reconnection, which was shown to be consistent with two-dimensional simulations [H, In this 
Letter I report three-dimensional high resolution simulations of the nonlinear evolution of the thin 
current sheet which spontaneously evolves into a turbulent current layer and shows a constant 
reconnection rate of ~ O.OISwa and the dissipation rate per unit area ~ O.OOGpw^, independent 
of the Lundquist number. Unlike the two-dimensional case, where reconnection is dominated by 
the ejection of large plasmoids and is very unsteady, I observed surprisingly steady rate. The 
turbulent current layer is slowly eating through the mostly undisturbed fiuid, with turbulence being 
fueled by the free energy of the oppositely directed magnetic fields. Around 40% of this energy is 
dissipated and 60% is converted into the turbulent fields. I conclude that the reconnection rate and 
the heating rate per unit area in the nearly-ideal current layers has a robust lower limit due to its 
inherent stochasticity. 



Magnetic X-points often collapse into thin current 
sheets [lj,ll5[- The Sweet-Parker (SP) model assuming 
resistive magnetohydrodynamics (MHD) predicts very 
low reconnection rate for most astrophysical and space 
magnetic configurations with high Lundquist number S. 
Moreover, the SP current sheet width is often below ion 
Larmor or ion skip depth di scales, which makes the 
model inapplicable in these situations. It was pointed 
out, however, that SP sheets are unstable above a cer- 
tain critical S* ~ 10^ [IS El; also, the plasma simulations 
reveal that moderately coUisionless systems are also un- 
stable The two-dimensional (2D) MHD results indi- 
cate that the reconnection speed is around 0.01 -^O.OSua 
12 , 13] and is due to hierarchical formation and ejection 
of plasmoids [l3| • Whether the same mechanism operates 
in the three-dimensional (3D) case was not clear. The 
reconnection due to ambient turbulence and the corre- 
sponding stochasticity of magnetic field lines [?[ gives the 
"fast" reconnection, independent on Lundquist number, 
but dependent on the level of ambient turbulence. Here 
I am investigating the high-^* case with current sheet be- 
ing unstable to tearing mode. In the absence of ambient 
turbulence the reconnection rate can also be fast as tur- 
bulence will be fueled by the reconnecting field. This 
could be understood as a lowest-stochasticity limit of the 
thin nearly-ideal current sheets. Magnetospheric obser- 
vations which show an enhanced level of turbulence inside 
current sheets [ll,[l3|, are supportive of this spontaneous 
reconnection picture. 

I simulated thin current sheets evolving into the turbu- 
lent current layers before the outflow becomes important. 



which is designated as the regime HI on Fig.[TJ Due to the 
fact that the outflow was unimportant, such a simulation 
could be realized in a very simple all-periodic setup, as 
shown on Fig. [2] I used simple pseudospectral code which 
has a number of advantageqj. Normal viscosity and mag- 
netic diffusivity rjV^B was used in one series of simula- 
tions and the hyper-diffusivity of the forth order, rj^V^B 
in the other. I defined the normal Lundquist number in 
terms of the box size L as S = Vj^L/rj and the hyper- 
Lundquist number as ^4 — vaL^/ti^. The Lundquist 
number S that gives the same ratio L/l^i as in the hyper- 
diffusive case could be expressed as 5 = S^^^^ , assuming 
Kolmogorov scaling. I used two dissipation schemes to 
check the influence of the dissipation to the bulk recon- 
nection and dissipation rates. As I will show below, the 
turbulence in the current layer seems to be local-in-scale, 
so at the sufficiently large Lundquist numbers the dy- 
namics of large scales, which determines the bulk recon- 
nection rate, is disconnected from the dynamics on small 
scales. I used unity magnetic Prandtl number Pr = ly/rj, 
so that the simulations could also be considered, roughly, 
large-eddy. The cubic resolutions was used to minimize 
the mutual inffuence of the two current sheets. 

The simulations were seeded with small initial pertur- 
bations, ~ 10^^ of the magnetic energy, and the bulk 
of the volume was almost undisturbed. The boundary 



^ A pseudospectral simulation can be "well-resolved" as long as the 
dimensionless maximum wavenur 
(ry^/e)^/* is the dissipation scale. 
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FIG. 1. A cartoon of high-Lundquist number magnetic re- 
connection. A generic magnetic x-point, (I), collapses into 
a thin current sheet, (II), which goes unstable and produces 
turbulent current layer (III), this layer expands with the re- 
connection speed Vr until it develops an outflow with speed 
Vout and the quasi-stationary state (IV) is reached with the 
inflow speed Vin = Vr is balanced by the outflow so that the 
current layer width A = Lvr/vout, neglecting compressibility. 



between the undisturbed volume, which had nearly zero 
current and the turbulent current layer was fairly well- 
defined, as shown on Fig. [21 The beginning of the current 
layer was determined as a point at which the current ex- 
ceeds a certain threshold in magnitude. The outflow was 
absent, so the reconnection rate was determined as the 
growth of the reconnection layer width, Vr = dA/dt. The 
measurements of the reconnection rate are presented on 
Fig. [3l The difference in the initial evolution of hyper- 
and normal-diffusion cases was due to faster tearing in- 
stability in the hyper case. I performed several simu- 
lations with different ratios Bzo/ByQ and different ByQ 
magnitudes. The reconnection rate was mainly propor- 
tional to VAy, while the dissipation to This is not 
surprising, since such dependence is suggested by units 
in the absence of dependence on S. The Vr dependence 
on S was studied by performing simulations with differ- 
ent resolutions, with S in the range of 1.7 6.4 • 10"* 
and 5'4^^° in the range of 2.5 5.4 • 10^. The evolu- 
tion of the current layer width and the inferred recon- 
nection rate are shown on Fig. [3] The three simulations 
2.5 • 10"' with different level of imposed 



with S^^^^ 

field ,620 = {0.5, 1.0, 2.0}i?j,o showed reconnection rates 
of {2.14, 2.10, 1.87}10~^VAy Such a weak dependence on 
BzQ is not surprising, as argued below. 

The energy density in code units, w = B^/2 + pv'^/2, 
consisted of the imposed field energy i?^Q and other con- 
tributions. The Bz flux through the current layer was 
conserved due to the conservation of the total B^ flux. 
The convenient dimensionless free energy of the recon- 
necting field could be defined as W ^ {2w — Bzo)/ByQ, 
which was unity in the undisturbed region. After the de- 
velopment of turbulence in the current layer, W ~ 0.6. 




FIG. 2. The all-periodic box reconnection setup with two 
current layers between the regions of reconnecting fleld ±Byo. 
The imposed mean field is B^o, the magnitude of B is shown 
on the surface of the box. 



This corresponds to the dissipation of 40% of the free en- 
ergy. The dissipation rate per unit area, therefore, was 
around 0.006pv\y. Most of the contribution to the re- 
maining 60% of W in the current layer came from B^ 
and By, while the contribution from SB^ = B^ — B^o 
was only around 2-4%. The contribution from kinetic 
energy was around 1% (normal diffusion) and 2% (hyper- 
diffusion). The turbulence in the current layer was 
anisotropic with respect to Bzo direction. This suggests 
that Bx and By components, carrying most of the en- 
ergy, represented Alfvenic perturbations, while the sub- 
dominant SBz was the slow-mode (pseudo-Alfven) per- 
turbation. It is, therefore, not surprising that the re- 
connection rate only weakly depends on Bzo, since the 
anisotropic Alfvenic turbulence is well-described by so- 
called reduced MHD, which possess the exact rescaling 
symmetry with respect to Bzo, see, e.g., This offers 
support to the models that investigated reconnection due 
to Alfvenic motions [7|, |2l| . It also explains why plasma 



simulations, such as |16j were similar to fluid simulations, 
as the reduced MHD is wel l-ap plicable to collisionless 
plasmas on scales above tl [23] ■ Despite the conserva- 
tion of flux, Bz was not constant in the current layer, so 
that Bz in the center was often as high as 1.3Bzo, while 
Bz on the edge was as low as Q.QBzo- This is probably 
due to the diamagnetism of turbulence, which is stronger 
on the edge, where turbulence is more intense, so that 
the diamagnetism has been pushing Bz flux towards the 
center. This interesting effect might have consequences 
for the particle acceleration in the current layers. 
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FIG. 3. The time evolution of tlie current layer width A and 
the inferred reconnection speed Vr (inset). The error bars of 
the width were obtained by varying the current thresholds so 
that the lower and upper thresholds differ by a factor of two. 



The power spectra E{^k1 + k^) of turbulent pertur- 
bations resembled so-called decaying MHD turbulence. 
The magnetic energy was dominant over kinetic energy 
on large scales, but the spectra converged with the small- 
scale magnetic and kinetic perturbations being almost 
in equipartition. The total energy spectral slope was 
around — 1.5-; — 1.7, roughly consistent with Kolmogorov 
or Goldreich-Sridhar [23j scaling. The spectral slopes be- 
tween -1 and -3 are indicative of local- in-scale turbulence, 
see, e.g., [3]. If and di are both smaller than the equi- 
librium layer width, 0.015L, the three-dimensional tur- 
bulence is local in scale and the bulk quantities, such as 
the average reconnection rate or the average dissipation 
rate could not be affected by plasma scales. However, the 
actual dissipation mechanism will, of course, depend on 
the properties of plasma. The difference between quiet 
dissipation in the solar wind that mostly heats plasma 
and the violent dissipation in the reconnecting current 
sheets in the solar corona, which produces electron ac- 
celeration, is the value of the dissipation rate per unit 
volume. When reconnection transitions between regimes 
II and III of Fig. [TJ the volumetric dissipation could be 
as high as e = Q.QQ&pv\/lsi>- This could push dissipa- 
tion scale Ir^ ~ g-i/4 below plasma scales and result in 
electron acceleration. 

It is interesting whether the measured reconnection 
rate will be affected by the presence of an outflow, as 
in the regime IV. When I increased the box size, the tur- 
bulent fluid motions had larger scales which could mimic 
the global outflow locally, but this did not increase the 
reconnection rate. Also, most of the activity which re- 
sults in a growth of the current layer happens on the 
boundary. Further 3D simulations with outflow will be 
able to test this, as soon as such simulations will break 



through the critical S barrier. 3D simulations of spon- 
taneous reconnection were performed in 25[ , but neither 
reconnection nor dissipation rates were reported. The 
minimum reconnection rate of 0.03w^ measured in Q is 
compatible with our results, as we did not drive ambient 
turbulence. 

The universality of the current layer large scale dy- 
namic properties is reminiscent of the properties of tur- 
bulent boundary layers in fluid dynamics. In the high-S* 
limit we expect all macroscopic properties to be expressed 
in terms of the p, va and L and independent of rj. This is 
similar to the other MHD phenomena, such as the non- 
linear small-scale dynamo 26|. 
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